Highly-efficient noise-assisted energy transport in classical oscillator systems 
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Photosynthesis is a biological process that involves the highly-efficient transport of energy cap- 
tured from the sun to a reaction center, where conversion into useful biochemical energy takes 
place. Even though one can always use a quantum perspective to describe any physical process, 
since everything follows the laws of Quantum Mechanics, is the use of quantum theory imperative 
to explain this high efficiency? Making use of the quantum-classical correspondence of electronic 
energy transfer recently introduced by Eisfeld and Briggs [Phys. Rev. E 85, 046118 (2012)], we 
show here that the highly-efficient noise-assisted energy transport described by Rebentrost et al. 
[New J. Phys. 11, 033003 (2009)], and Plenio and Huelga [New J. Phys. 10, 113019 (2008)], as 
the result of the interplay between the quantum coherent evolution of the photosynthetic system 
and noise introduced by its surrounding environment, it can be found as well in purely classical 
systems. The wider scope of applicability of the enhancement of energy transfer assisted by noise 
might open new ways for developing new technologies aimed at enhancing the efficiency of a myriad 
of energy transfer systems, from information channels in micro-electronic circuits to long-distance 
high-voltage electrical lines. 

PACS numbers: 87.10.-e, 82.20.Nk, 82.20. Rp 
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Because of its undoubted importance for all life on 
earth, molecular mechanisms of energy transport in pho- 
tosynthetic light-harvcstingcomplexes have been a sub- 
ject of study for decades [l|-|4j. In recent years, a re- 
newed interest on this topic has arisen Q , mainly due 
to the unexpected observation of long-lived electronic co- 
herences in the energy transfer process of photosynthetic 
systems, particularly in the nowadays most widely inves- 
tigated system, the Fenna-Matthcws-Olson (FMO) com- 
plex 0-Q. 

As a consequence of these findings, several theoretical 
studies have been devoted to describing how coherence 
effects in a quantum scenario might play an important 
role in the remarkably high efficiency of energy transfer in 
photosynthetic systems [ld - IHa ] . This is specially notable 
since it takes place in a scenario apparently not propitious 
for the observation of quantum effects. In particular, it 
has been suggested that high efficiency transport arises as 
a result of the dynamical interplay between the quantum 
coherent evolution of the photosynthetic system, and the 
dephasing noise introduced by its surrounding environ- 
ment, a phenomenon called environment-assisted quan- 
tum transport (ENAQT) fill ] or dephasing-assisted en- 
ergy transport [141 ]. 

As stated in Ref. ENAQT can be understood 

as the suppression of coherent quantum localization me- 
diated by noise, helping the excitation to move faster 
through the photosynthetic system, thus increasing the 
efficiency of energy transport. In this way, ENAQT 
might be seen as a phenomenon that exists only in a 
regime where the quantum and classical worlds overlap. 
Notwithstanding, making use of the quantum-classical 
correspondence of electronic energy transfer presented in 
Ref. lift, we show here that the same effect can also be 



found in purely classical systems. Our departure point 
is based on the consideration that aggregates of coupled 
monomers (such as the FMO complex) can also be de- 
scribed as a system of weakly-interacting classical oscil- 
lators [l7j . We then demonstrate that the noise-assisted 
enhancement of transport efficiency in the FMO complex, 



shown in Refs. [I3|,ll4| , and based on a pure quantum for- 



malism, it can also be found in a purely classical model, 
without the need to resorting to quantum effects. This 
result arises as a consequence of the gen eral quantum- 
classical equivalence described in Ref. [161 ] . 

For the sake of comparison and clarity, we will first 
model the FMO complex as a quantum system of N in- 
teracting sites, where the interaction of each site with its 
surrounding environment is modeled by a pure dephas- 
ing process. We have adopted this model because of its 
extended use for describing noise-assisted energy trans- 
fer processes in photosynthetic systems [H El. Next, 



we will present the classical model of Refs. [16|, |17| . 
which corresponds to a system of N weakly-coupled har- 
monic classical oscillators. In this case, environmental 
effects are introduced by assuming that the frequency of 
each oscillator vary stochastically as a Gaussian Markov 
process. Finally, we will solve both models using the 
site energies and coupling coefficients for the FMO com- 
plex of Prosthecochloris aestuarii to show that the same 
environment-assisted energy transfer effect can be found 
in both, the classical and quantum models. 

The Hamiltonian of a system comprising N interacting 
sites in the presence of a single excitation is given by 
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where |n) denotes the excitation being at site n. The 
nth-site energies and the coupling between sites n and m 
are described by e n and V nm , respectively. 

We make use of a simple model where the dynamics of 
the system interacting with a surrounding environment 
is described by a Lindblad master equation, which in the 
Born-Markov and secular approximations writes as (l8j 
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dt 
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Here, the interaction of the system with the envi- 
ronment is characterized by a pure dephasing pro- 
cess given by the Lindblad operator Cdeph [p\ nm = 
- [1/2 (7„ + 7,„) - y/Jn'jmSnm] p nm , with 7„ being the 
dephasing rates. Although the pure dephasing model is 
not able to capture important aspects of electronic en- 
ergy transfer, such as phonon relaxation it provides 
a useful description of environmental effects in a simple 
way. To quantify the transfer of energy from a chosen 
site k to the reaction center, we have phenomenologi- 
cally introduced an irreversible decay process (with rate 
T) described by the operator D, which is given by [2^ 
D[p\nm = - r {l fc ) ( fc UL m > wher e {...} stands for the 
anticommutator. 

Making use of equation ([2]), one can define a measure 
for the efficiency of energy transport as the population 
transferred to the reaction center, within a time t, as 



e // = 2T / (k\p(s)\k)ds. 



(3) 



Equations (|2|) and ((3J constitute the quantum equations, 
which have to be compared with the equations that will 
be obtained in the classical model. 

For the classical case, we consider an ensemble of N 
coupled harmonic oscillators, each with mass M and fre- 
quency u) n . The temporal evolution of the system is de- 
scribed by a classical Hamiltonian, which in terms of the 
position q n and momentum p n of each oscillator, reads 
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where K nm stands for the coupling coefficient between 
the oscillators. By defining a new dimcnsionless com- 
El: 



plex amplitude 



z n (t) = q n (t) + ip n (t), with q n 



(Muj n /2hyl 2 q n and p n = (2hMuj n ) 



-1/2 



Pn • 



the Hamil- 



ton equations of motion of the system can be cast into a 
single equation 



dz n 
dt 



(5) 



Re{...} stands for the real part of a complex number, 



and K,„ 



K nm /(My/UJ n UJ m ). 
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FIG. 1. Arrangement of the BChla molecules of a single unit 
of the Fenna-Matthews-Olson (FMO) complex. The figure 
was created using PyMOL [3^], and is based on the PDB 
entry 3ENI. 



To include environmental effects, we proceed in the 
same manner as in the construction of a Kubo oscilla- 
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For this, we assume that the frequency of 
each classical oscillator varies randomly as a stochastic 
process: ia n (t) = u> n + 4> n (t). uj n is now the average fre- 
quency of the nth oscillator and <fi n (t) describes a Gaus- 
sian Markov process with zero average (Wiener process), 
i.e., (</>„(£)) = and (4> n (t) cj) m (t')) = 7„<5 nm <S (t - t'), 
where (...) denotes stochastic averaging. 

In Ref. , it has been shown that one can transform 
Eq. (|5|) , within the framework of Ito calculus 2J] , into 
a classical master equation that describes the temporal 
dynamics of the system of coupled harmonic oscillators, 
when interaction with the surrounding environment is 
taken into account. To describe the transfer of excitation 
from the fcth oscillator to the reaction center, we extend 
this result, and introduce an irreversible decay process 
(with rate T), described by V [a] nm =-T{\k)(k\, a} nm , 
where a nm = (z n z^). In this way, we can write the 
classical master equation as 



da r . 



dt 



= n[a] nm + £[a] nm +V[a} 7 
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with V nm = K nm h/2, and 

— T {VnjG'jm ~ VjmCnj) j 
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function of the dephasing rate 7 obtained by means of 
the quantum equations ([2]) and ([3]). Notice that at low 
dephasing, i.e., with environment effects not considered, 
coherent evolution of the system leads to an efficiency of 
about 90%. When increasing the dephasing, efficiency 
grows to almost 100%, showing that the environment af- 
fects the system in such a way that it becomes more 
efficient for transferring energy to the reaction center. 
Finally, for stronger dephasing, efficiency drops rapidly 
and almost no energy is transferred to the reaction cen- 
ter. Qualitatively similar results have also been obtained 
for the case of the FMO complex of Chlorobium tepidum 



FIG. 2. Energy transfer efficiency as a function of the dephas- 
ing rate 7 obtained from the quantum-mechanical equations. 
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The energy transfer efficiency within the ensemble of os- 
cillators is then given by 
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a (t) = a (t) I ^2 n a nn is the normalized classical density 
operator. 

Equations © and © represent the classical equations, 
whose results have to be compared with their quantum- 
mechanical counterpart, equations ([2]) and ([3]). To this 
end, we can make use of the site energies and coupling 
coefficients for the FMO complex of P. aestuarii [25| . 
The FMO is a pigment-protein complex that guides the 
energy from the light-harvesting chlorosomes to the re- 
action center in green sulfur bacteria [2(| 27 1. It is 



a trimer of three identical subunits interacting weakly 
with each other. Each subunit is composed of seven 
bacteriochlorophyll-a (BChla) molecules embedded in a 
scaffolding of protein molecules, as shown in Fig. 1. The 
FMO complex is generally modeled by a network of seven 
different sites, where the dynamics of a single excitation 
through the complex is governed by the specific values of 
the site energies (e n ) and the coupling coefficients (V nm ). 
In particular, we will use the values of the site energies 
and coupling coefficients for P. aestuarii, as stated in 
Tables 2 and 4 of reference [25| . The initial state of the 
system corresponds to a single excitation in site 1. In 
the FMO, the BChl 3 is in the vicinity of the reaction 
center [25|. Thus, we take this site (k = 3) as the main 
excitation donor to the reaction center, with a transfer 
rate estimated to be Y = 1 ps -1 (28j . Furthermore, for 
the pure dephasing process, we consider that dephasing 
rates are the same for all sites (7 = 7«) and that the effi- 
ciency of energy transfer is limited by the finite excitation 
lifetime (t ~ 1 ns). 

Figure 2 shows the efficiency of energy transfer as a 



We now turn our attention to the case of the classical 
model by solving equations ^ and ©• Figure 3 shows 
the efficiency of the energy transfer as a function of the 
dephasing rate. We observe that the same noise-assisted 
effect is also present in the purely classical model. For 
the sake of comparison, Fig. 3 also shows the solution 
of the quantum mechanical model (dashed line). Notice 
that both solutions agree for dephasing rates up to 10 3 
ps -1 . However, for larger values of dephasing the quan- 
tum and classical solutions differ from each other. This 
is in agreement with the fact that both solutions are the 
same, provided that the condition 7 <C uj n is satisfied 

Noise-assisted energy transport in disordered systems 
has been understood as the suppression of coherent quan- 
tum localization through noise, bringing the detuned 
quantum levels into resonance and thus facilitating the 
energy transfer [3, EH. Notwithstanding, the results 
presented here show that the same effect can also be 
found in purely classical systems. This implies that one 
can make use of such systems in order to simulate the 
intricate energy transfer mechanisms that take place in 
molecular aggregates, such as the photosynthetic FMO 
complex. 

Recently, it has been suggested that classical LC cir- 
cuit oscillators (where L stands for inductance and C 
for capacitance) can be used to model coupled quantum 
two- level systems (2^|. Hence, one could devise an ex- 
perimental apparatus comprising eight electrical oscilla- 
tors with the eighth acting as the reaction center, which 
would be strongly coupled to one of the remaining oscil- 
lators. Then, by stochastically modulating the frequen- 
cies w„, and properly controlling the noise intensity 7„, 
one would be able to observe the noise-assisted energy 
transfer phenomenon by monitoring the signal present in 
the eighth oscillator. These classical simulations could 
be further used to compare with the recent experimen- 
tal proposal of noise-assisted transport based on coupled 
quantum-optical cavities [30j ] . 

The concept of noise-assisted energy transport has 
been extensively used for describ ing the inner working 
of quantum and classical systems [311 ]. Along these lines, 
the particular enhancement effect described in this pa- 
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FIG. 3. Energy transfer efficiency as a function of the dephas- 
ing rate 7 obtained from the classical equations (solid line). 
For the sake of comparison, we have included the results for 
the quantum equations (dashed line). 



per might open a new research direction towards new 
methods for enhancing the efficiency of a myriad of en- 
ergy transport systems that inevitable live in a noisy 
environment, from small-scale information and energy 
transfer systems in microwave and photonic circuits, to 
long-distance high-voltage electrical lines. In this way, a 
specific feature initially conceived in a quantum scenario 
(environment- assisted energy transport) is shown to arise 
as well in a purely classical context, widening thus the 
scope of possible quantum-inspired technological applica- 
tions. 

To conclude, the search and demonstration of systems 
where to observe quantum-mechanical effects with no 
classical counterpart is a subject of lively interest and 
debate (32-34|. Biological systems arc not, in principle, 
a propitious scenario for the observation of quantum fea- 
tures, such as quantum superposition, interference or en- 
tanglement. Nevertheless, one can always use a quantum 
perspective to describe any physical process, since ev- 
erything follows the laws of Quantum Mechanics. This 
does not mean, however, that in certain cases a purely 
classical model may not similarly reproduce some of the 
results predicted by the full quantum-mechanical model, 
since classical physics emerge, after all, from quantum 
physics under many circumstances. 
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